Preprint APCTP/98-22 
KEK-TH-594 



Rep-ph/ 9810369 

Decay Constants and fr> s 

from B° -> D+Z-p and B° -+"d + DjW Decays 

Yong-Yeon Keum [] 

Asia Pacific Center for Theoretical Physics 
207-43 Cheongrangri Dong, Dongdaemun-Gu, Seoul 130-012, Korea 

Theory Group, KEK, Tsukuba, Ibaraki 305-0801 Japan 

Abstract 

We calculate the decay constant of D s and D* with B° —> D + £~v and B° — > D+DJ" 
decays. In our analysis we assume the factorization ansatz and use two different form 
factor behaviours (constant and monopole-type) for i*o((? 2 ). We also consider the QCD- 
penguin contributions in hadronic decays within the NDR renormalization scheme in a 
NLO calculation. 

We estimate the decay constant of the D s meson to be 219 ±46 MeV for (pole/pole)- 
type form factor and 239 ± 50 MeV for (pole/constant)-type form factor. For D* meson, 
we predict /d* = 336 ± 79 MeV, and /d*//d s = 1-41 ± 0.41 for (pole/constant)-type 
form factor. 
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1. Introduction 

Measuring purely leptonic decays of heavy mesons provides the cleanest way to determine 
the weak decay constants of heavy mesons, which connect the measured quantities, such as the 
BB mixing ratio, to CKM matrix elements V c b, V u b- However, at present it is not possible to 
determine /b,/b s fo s and fp* experimentally from leptonic B and D s decays. For instance, 
the decay rate for is given by 

T(Dj - fu) = ^f Ds m]M Ds (l - j±J \V cs f (1) 

Because of helicity suppression, the electron mode —>■ e + u has a very small rate. The relative 
widths are 10 : 1 : 2 x 10 for r + z/, fi + u and e + u final states, respectively. Unfortunately the 
mode with the largest branching fraction, t + i>, has at least two neutrinos in the final state and 
is difficult to detect experimentally. So, theoretical calculations for decay constant have to be 
used. The factorization ansatz for nonleptonic decay modes provides us a good approximate 
method to obtain nonperturbative quantities, i.e., form factors and decay constants which are 
hardly accessible in any other way || [| . 

There are many ways that the quarks produced in a nonleptonic weak decay can arrange 
themselves into hadrons. The final state is linked to the initial state by complicated trees of 
gluon and quark interactions, pair production, and loops. These make the theoretical descrip- 
tion of nonleptonic decays difficult. However, since the products of a B meson decay are quite 
energetic, it is possible that the complicated QCD interactions are less important and that the 
two quark pairs of the currents in the weak Hamiltonian, group individually into the final state 
mesons without further exchanges of gluons. The color-transparency argument suggests that a 
quark-antiquark pair remains a state of small size with a correspondingly small chromomagnetic 
moment until it is far from the other decay products. 

Color transparency is the basis for the factorization hypothesis, in which amplitudes fac- 
torize into products of two current matrix elements. This ansatz is widely used in heavy quark 
physics, as it is almost the only way to treat hadronic decays. 

In this paper we consider a way to determine weak decay constants fo a and fo* under 
factorization ansatz. In section 2, we discuss how to extract the unknown parameter | V c bFf D (0) | 



1 



from the branching ratio of the semileptonic decay B° — > D + £u. In order to check the validity 
of the factorization assumption, we study the nonleptonic two body decays, for instance, B — > 
Dp, Dn and DK^ in section 3. In section 4, we discuss the determination of fjj a and fu* 
from B° — > D + D~^ decay modes. In our analysis, we consider the QCD-penguin effects which 
amount to about 11 % for B — ► DD S and 5 % for B — > DD*, which have been mentioned in 
Ref.l. 



2. Semileptonic Decay B° — > D + l ~v 

From Lorentz invariance one finds the decomposition of the hadronic matrix element in 
terms of hadronic form factors: 



<D + {p D )\J^\B\p B )> 



™B - m D 



BDf „2\ 



m u — m 



(2) 

where = c^^b and = (j>b — Pd)h- I n the rest frame of the decay products, Fi(q 2 ) and 
F Q (q 2 ) correspond to 1~ and + exchanges, respectively. At q 2 = we have the constraint 



(3) 



since the hadronic matrix element in (0) is nonsingular at this kinematic point. 

The q 2 distribution in the semileptonic decay B° — > D + l~v is written in terms of the 
hadronic form factor Ff D (q 2 ) as 



dg 2 24tt 3 
where the g 2 dependent momentum K(q 2 ) is given by 

1 



(4) 



K(q 2 



2m b 

2 



m B + m D — q ) — Am B m D 



1/2 



(5) 



In the zero lepton-mass limit, < q < (m# — % 
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For the q 2 dependence of the form factors, Wirbel et al. || assumed a simple pole formula 
for both Fi(q 2 ) and F (q 2 ) (we designate this scenario 'pole/pole'): 



F 1 {q 2 )= Fx(0)/(1 



m 



2 n 



F (q 2 )=F (0)/(l-^) 



m 



2 /' 

Fo 



(6) 



with the pole masses 



m Fl = 6.34 GeV, m Fo = 6.80 GeV. (7) 

Korner and Schuler |§ have also adopted the same q 2 dependence of Fi(q 2 ) and F (q 2 ) given 
by @ and (]?]). On the other hand, the heavy quark effective theory (HQET) gives in the 
mb )C — > oo limit the relation between Fi(q 2 ) and F (q 2 ) given by [11|, jl2| 

,2 



1 - 



q 



F 1 (q 2 



(8) 



(m B + m D ) 2 

The combination of (§) and @ suggests that -Fo(<? 2 ) is approximately constant when we keep 
the simple pole dependence for Fi(q 2 ). Therefore, in this paper, as well as the above 'pole/pole' 
form factors, we will also consider the following ones (designated 'pole/const.'): 

,2 



F 1 tf) = F 1 (0)/(l 



q 



m 



2 /' 
Fi 



F (q 2 ) = F (0), 



(9) 



with 



m Fl = 6.34 GeV. (10) 

By introducing the variable x = q 2 /m 2 3 , which has the range of < x < (1 — ^) 2 in the 
zero lepton mass limit, (H) is written as 

a 3 [i, % A 



dT(B i 



D+l-v) 



dx 



A[l, — x \ 



G F m% 

1927T 3 



V cb Fi 



BD 



m 



B 



:i + 



m 



D 



x 





m| 








1/2 


m 2 D 







X 



2 ■ 



(11) 



Then the branching ratio B(B° — > D + l v) is given by 



G F m 2 B 2 trb 2 



y/2 ' T B 192vr 2 
2.221 x 10 2 \V cb Ff D 



Kb ^i (0)| x I 



x L 



(12) 
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where the dimensionless integral I is given by 



m B 1 J 

ax 



4^1 



3/2 



1 



0.121 



(13) 



In obtaining the numerical values in (0) and (0), we have used the following experimetal 
results H: m D = m D + = 1.869 GeV, m B = m B o = 5.279 GeV, T B = T B o = 4.219 x 10~ 13 
GeV (t b o = (1.56 ± 0.06) x 10~ 12 s), and G F = 1.166 39(2) x 10~ 5 GeV -2 . Since B(B° -> 
D + rv) = (1.78 ±0.20 ±0.24) x 10~ 2 is obtained experimentally gj|, the value of \V cb F? D (0)\ 
can be extracted from ([T2]) . Following this procedure, we obtain 



V cb Ff D (0) \ = (2.57 ±0.14 ±0.17) x 10 



-2 



(14) 



In the calculations that follow, we will use \V c b Ff (0) \ = (2.57 ± 0.22) x 10 which combines 
the statistical and systematic errors in ([14]). 



3. Test of Factorization with B° — > D + p and B° — > D + tt , and Prediction 
of Branching Ratio B{B° -> D + K~M) 

We start by recalling the relevant effective weak Hamiltonian: 

H eS = ^WftPi(AO0i + C 2 (fi)0 2 ] + E.G., (15) 
O x = (dT p u)(cT p b), 2 = {cT p u){dT p b), (16) 

where Gp is the Fermi coupling constant, V cb and V u d are corresponding Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements and T p = 7 p (l — 75). The Wilson coefficients Ci(fi) and 
C 2 (fi) incorporate the short-distance effects arising from the scaling of 7i e R from \x = mw to 
fi = 0{m b ). By using the Fierz transformation under which V — A currents remain V — A 
currents, we get the following equivalent forms: 

C x Ox + C 2 2 = (C 1 + ^-C 2 )0 1 + C 2 {dT' , T a u){cT p T a b) 

= {C 2 + ^C l )0 2 + C l {cr p T a u)(dT p T a b), (17) 
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where N c = 3 is the number of colors and T a s are SU(3) color generators. The second terms 
in dl7|) involve color-octet currents. In the factorization assumption, these terms are neglected 
and TC e s is rewritten in terms of "factorized hadron operators" || llOfl : 

Gi 



H eS = ^V cb V: d (a 1 [^ p u} H [cT p b} H + a 2 [cr p u} H [^ p b} H ) + H.C. 



;is) 



where the subscript H stands for hadronic implying that the Dirac bilinears inside the brackets 
be treated as interpolating fields for the mesons and no further Fierz-reordering need be done. 
The phenomenological parameters ai and a 2 are related to C\ and C 2 by 



O-l — C\ + — C 2 , CL 2 — C 2 + J^~^l- 



(19) 



From the analyses of A.J. Buras |]I4[ , the parameters a\ and a 2 are determined in the Next- 
to-Leading-Order (NLO) calculation in the Naive Dimensional Regularization (NDR) scheme 

as 

ai = 1.02 ±0.01, a 2 = 0.20 ±0.05. (20) 
For the two body decay, in the rest frame of initial meson the differential decay rate is given 

by 



dr 



32tt 2 ' 1 M 2 



IpiI 



[(M 2 - (mi + m 2 ) 2 )(M 2 - (mi - m^ 2 )] 1 / 2 



2M 



(21) 
(22) 



where M is the mass of initial meson, and mi (m 2 ) and pi are the mass and momentum of one 
of final mesons. By using (0), (|l"8D and < 0\T tl \p(q,e) >= e^(q)m p f p , (|21|) gives the following 
formula for the branching ratio of B° — > D + p~: 

2 , f 2 

V y/2 J l ' uai 16tt T B ^mf 



m B 



-) 2 )(i-(- 



m B 



3/2 



13.25 x \V cb F* D (m 2 p )\ 2 . 



(23) 



In obtaining the numerical values in (|23|), we have used m p = m p + = 766.9 MeV, f p = f p + = 216 
MeV, V u d = 0.9751 [031, and the parameters given below (O). 
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For at we use the value given in fl2"0"|). Then, by using the formula ( |2"3"D with the values of 
\Vd>Ff D (Q)\ 2 , (F^ D (0) = Ff D (0)) given in (0, we obtain the branching ratio B(B° -> £> + p~) 
presented in Table 3. 

For the process _B° — > D + K*~, by using < 0|r At |i^*(g, e) >= e^(q)m K * f K *, we have 

3/2 



16vr r B 



m 



,m D + m^* 2 



) 2 )(i-( 



BD/^2 \|2 



- * 2 



0.67 x IK^i (mi 



(24) 



where we have used = m K - = 891.59 MeV, /#* = = 218 MeV, and V us = 0.2215 
|33|. By using (g|) with \V cb Ff D (0)\ 2 in (0), we obtain the branching ratio B(B° -> D+K*~) 



presented in Table [I]. 

By using (|2|), ( [IB] ) and < 0|r M |7r(g) >= iq^fx, ( pi]) gives the following formula for the 
branching ratio for the process B° — > D + tt~: 



B(B° -> D+tt" 



\ V2 ) 



2 <-2 

7T 

16tt T B m| 



IT/ 12 1 m B „2 /tt 



x 1 



,m D + m w 2 



m B 



) 2 )(i-( 



1/2 



5.42 x |K 6 ^K)| 2 . 



(25) 



where we have used m,,- = m,,-- = 139.57 MeV and = /„■- = 131.74 MeV fl33] . By using 
the formula (0) with the values of \V cb F^ D (0)\ 2 , (F^ D (0) = Ff D (0)) in ([14]), we obtain the 



branching ratio for B° 
For the process B° 

B(B° -> D+i<r 



_D + 7r presented in Table [l|. 

D + K~ , by using < 0|r M |K _ (g) >= iq^fx-, we have 

/ G F m 2 B \ 2 2 1 m B 2 pBfl^2 M2 



v V2 y 



x 1- 



m 2 D \2 



rn 



B 



16tt T B 1 m| 
m D + m K 2 



m B 



cbr \m K) \ 

m D - m K 2 



1/2 



o.4i x iK b F »)r, 



(26) 



where we have used = m^- = 493.68 MeV, = f K + = IQO.QMeV p3| . By using ( Effi) 
with | i 7 ^^ ( 0) | 2 in ([T4]) , we obtain the branching ratio B(B° — > D + K~) presented in Table |l[ 



An inspection of Table 1 shows that the factorization method works well in B° — > D + n , D + p 
decays. We predict branching ratios : 

B{B° -> L>+*T) ~ 2.7- 10~ 4 

B(B° -> D + K*-) ~ 4.6 • 10~ 4 (27) 
which is certainly reachable in near future. 



4. Determination of /p* and f D from 5° — ► D + D * and 5° — > D + D~ 



(a) At Tree level 

In this section we calculate the vector and pseudoscalar decay constants fo* and fo 3 from 
the experimentally obtained branching ratios of the the exclusive decays B° — > D + D~* and 
50 ^ D+£) - By using (ii) and < |r^|D;( g ,e) >= e^q)m m f DS , gives the following 
formula for the branching ratio of _B° — > D + D~*: 



B(B° -> D+D; 



(Gprn^Y 2 2 1 m B 2 f 2 

,m D + m D * 2 



•'"-d „2 J-D; I rpBDf 2 \|2 



r i-( 



1.32 • 10 2 GeV 2 ) • / 2 , • |K fe Ff D (m^)| 2 



3/2 



(28) 



where we have used mo* = m D -* = 2112.4 MeV V cs = 0.9743 [33|. Then, from (|28D we get 



/d* = (0.87 x lO" 1 GeV) 



B(B° -> 



(29) 



For the process S° — ► D + D s , by using < 0|r At |D s (g) >= iq^fDs, we have 

sf-B° — ► d + dt) = f^#Vi^ri^rT^^ a i^i F o BD (<)i 2 

1 7,1 ' 167r 1 b m B 



x(l 



m 2 B ' 



,m D + m fl 2 



) 2 )(i-( 



m D - m D , 



s \2 



m B 



1/2 



(2.45 -10 2 GeV 2 )-f Ds -\V cb F BD (mlJ\ 2 , 



(30) 
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where we used mn 



m 



D- 



1968.5 MeV 



. From 



we get 



f Ds = (0.64 x GeV) ■ 



B{B° -> D+D; 



\V cb F BD (mlJ\ ■ 

Browder et al. [15] present the following experimental results for the branching ratios: 



(31) 



B(B° -> 

-6(5° -> £> + £>7) 



'1.14 ±0.42 ±0.28) x 10" 



;i.l4±0.50) x 10" 



(0.74 ± 0.22 ± 0.18) x 10~ 2 = (0.74 ± 0.28) x 10" 



(32) 



where we have combined the statistical and systematic errors. By using ([32] ) and the values 
of \V cb Ff D {m 2 D ,)\ and \V cb F BD (m 2 Ds )\ given in Table 2, from (0) and © we obtain the the 
following results: 



f D * = 322 ± 76 MeV, f Ds = 196 ± 41 MeV for (pole/pole), 
fjj* = 322 ± 76 MeV, f Ds = 214 ± 45 MeV for (pole/const.) 

The ratio of the vector and pseudoscalar decay constants fjj* / fjj s is given by 



(33) 



f. 



1.36 x 



\V cb F BD (ml s 
\V A F? D (ml. 



x 



B(B Q 



D + D: 



B(B° -> D+D; 



1/2 



(34) 



which gives 



D: 



fD s 



1.64 ±0.48 for (pole/pole), 
1.51 ± 0.44 for (pole/const.). 



(35) 



(b) Including Penguin contribution 

The effective Hamiltonian for AB = 1 transitions is given by 

H eS = ^tKaVyCOT + C 2 O u 2 ) + V cb V: q [C x O\ + C 2 O c 2 ) - V tb V* £ CA], (36) 

V2 i=3 

where q = d,s and Cj are the Wilson coefficients evaluated at the renormalization scale /i, and 
the current- current operators 0\ ,c and 0%' c are 

O" = {u a b a ) v _ A (q f) up) v _ A 0\ = {c a b a ) v _ A (qpCp) v ^ A 

0% = (u{3b a ) V - A (q a up) V - A 0% = {cpb a )v~A{q a Cf3)v~A, (37) 



and the QCD penguin operators O3 — 0§ are 



3 = {q a b a )v-A Y.Wp)v-A 4 = {qpb a )v_ A ^{qW^v-A 

q> q> 

5 = {.q a b a )v-A^2%q'a)v+A 6 = 



a)V-A ' V 



(38) 



In ( |36l) we have neglected the effects of the electroweak penguin operators and the dipole 
operators, since their contributions are not important for the calculations of this paper. 

When we take m t = 174 GeV, m b = 5.0 GeV, a s {M z ) = 0.118 and a em {M z ) = 1/128, the 
numerical values of the renormalization scheme independent Wilson coefficients Ci at \l — m b 
are given by [[L6[] 



Ci = -0.3125, C 2 = 1.1502, 

C 3 = 0.0174, C A = -0.0373, C 5 = 0.0104, C 6 = -0.0459. 
The effective Hamiltonian in ([36]) for the decays B° — > D + D^^ can be rewritten as 

h cS = %\v cb v: s (cfoi + cfo c 2 ) - v tb v t ; f; cf oj, 

V2 i=3 

where Cf ff are given by |I7| 



(39) 



(40) 



with 



Ceff /~feff /m /^/cfT ID / AT /^eff /~i i D 

1 — tvl? ^2 — ^2) ^3 — ^3 — -n s/iV c , o 4 — O4 + r s . 

Cf = C 5 -P„/N c , Cf = C 6 + P B 



P s = ^-G(m q , q 2 ,»))CM, 

2 f 1 — x(l — x)g 2 

G(m qj q , a) — —4 / a;(l — a;) ln( — = )d:r, 

jo /r 



(41) 



(42) 



where g denotes the momentum of the virtual gluons appearing in the QCD time-like matrix 
elements, and iV c is the number of colors. Assuming q 2 = ml/2, we obtain the analytic formular 
for G(m q , q 2 , fi) : 

,2 



GK, -, /,) = --In + 3 



in 



1- VT=|/ 



1 + vT^ 



+ 27T 



(43) 



with y = 9>m? q /m 2 . By inserting the values for m q = m c = 1.45 GeV, and N c = 3 into ( f4"T| ) and 
(f|3|), we get the values C e ^ [i — 1 ~ 6) which is shown in Table |2|. 

The decay amplitude ^4t+p(-B° - ► D + D~) =< D + D~\H e s\B° > is given as follows: 



A T+P (B° ^ D + D;) = ^L[V cb V c * s a 2 -V tb V t :(a 4 + 2a 



ni 



V2 

Gp 
V2 



'(m b -m c )(m s + m c ) 



)] M a 



V cb V*a 2 



m 



1 + ^1 + 2^ 

a 2 a 2 (m b - m c ) (m s + m c 



Mr. 



where 

M a =< D;\s^ l5 c\0 >< D + \c %l b\B° >= -if Ds (ml - m 2 D )F BD (m 2 Da 
On the other hand, we have 

Gp . 



(44) 



(45) 



A T+P (B° -> D + D s 



V2 



[V cb V*a 2 - V tb V*a 4 ] Mi 



^V cb V; s a 2 (l + ^) M b 



(46) 



where 



M b =< D* s \sYl 5 c\0 >< D + \c^b\B° >= m Dt f Dt [e(q) ■ (P B + P D )]F( 



nr 



D: 



(47) 



We can estimate the penguin contributions for each process from ([0]) and ( flTf ) as follows : 



For B° -> D + D: 



For B 



o 



D+D*~; 



A P 




0,4 


At 




a 2 


A P 




04 


A T 




a 2 



rn 



Ds 



a 2 (m b - m c ){m c + m s 
4.5% 



10.9% (48) 
(49) 



Here we have used the values for m c = 1.45 GeV, and m s = 170 MeV. The penguin contribu- 
tions for B — > DD S is two times larger than that in B — > DD*, and is a sizable effect. This 
results agrees well with that of Ref. |4j]. Hence penguin contributions can affect the value of the 
decay constants fo a and fp>*. 

From fl3~T|), (|2"9"|), ([P]) and (^) we obtain the following results including penguin contribu- 
tions : 



Id* — 336 ± 79 MeV, fp> s = 219 ± 46 MeV for (pole/pole), 
f D * = 336 ± 79 MeV, f Ds = 239 ± 50 MeV for (pole/const/ 
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(50) 



From (0), ( 
given by 

which gives 



and ( |4"6j ) the ratio of the vector and pseudoscalar decay constants /d*//d s is 



D: 



1.36 



D a 



\V cb F BD (ml s )\ 



£(5 C 



D + D~*) 



1/2 



0.8951 
0.9582 



(51) 



/d; 



1.53 ±0.45 for (pole/pole), 
1.41 ± 0.41 for (pole/const. 



(52) 



The decay constant depends the q 2 behaviour of the form factor Fo(q 2 ). However the amount 
of change is less than 10% as shown in (0). This shows that the decay constant is not strongly 
dependent on the behaviour of the form factor. As shown in Table ^], when we include the 
QCD-penguin contributions, the value of the decay constant fu* is increased by 4%; however, 
for fu s it is increased by up to 12%, and the ratio /dj//d s is decreased by 6%. Table 3 shows 
that our result for fn a agrees well with others' predictions and experimental data, except that 

For the ratio /d*/ /d s , our results have a value greater than 1. On 
and Hwang and Kim 2 [TOl have a value less than 1, which is 



of Hwang and Kim 



the contrary, Browder et al. 
quite strange, because their results show that the decay constant of the vector meson is smaller 
than that of the psudoscalar meson with the same quark contents. In fact, the decay constant 
of p meson is 1.5 times greater than that of ir meson. By measuring the ratio fn*/fD s , we can 
check directly their analysis and methods. 



5. Conclusion 

Within the factorization approximation, we have calculated the weak decay constants fo s 
and fo* from the semileptonic decay, B° — > D + £~v and the hadronic decay modes, B° — > 
D + D < f\ We have also considered the effect of two different g 2 -dependences of the form factor 
F^ D (q 2 ). The value of fu a is altered by less than 10% for different form factors. In our analysis, 
we also considered the QCD-penguin contributions in hadronic two body decays within the 
NDR renormalization scheme in a next-to-leading order calculation. The penguin effects for 
B — > DD S decay are quite sizable. We obtain fo a = 219 ± 46 MeV for the monopole type of 
F BD , f Da = 239 ± 50 MeV for the constant F BD . We obtain f D , = 336 ± 79 MeV for the 
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D s * meson. These predicted values will be improved vastly when the large accumulated data 
samples is available at Belle and BaBar experiment in near future. 
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B(B° -> D+p-) 


B(B° -> D+K-*) 


B(B° -> D+7T-) 


-6(5° -> T>+^-) 




xlO 3 


xlO 4 


xlO 3 


xlO 4 


(pole/pole) 


9.01 ± 1.54 


4.62 ± 0.79 


3.58 ±0.61 


2.74 ±0.47 


(pole/const.) 


9.01 ± 1.54 


4.62 ±0.79 


3.57 ±0.61 


2.71 ±0.46 


Experiments 


8.4 ± 1.6 ±0.7 




3.1 ±0.4 ±0.2 





Table 1: The obtained values of the branching ratios and existing experimental values. 



Coefficients 


Real Part 


Imaginary Part 




- 0.313 


0.0 


cs" 


1.150 


0.0 


cf f 


2.20 T0~ 2 


5.13 TO- 3 


cf f 


-5.11 T0~ 2 


-1.54 TO -2 


ci" 


1.50 TO" 2 


5.13 TO -3 


eg" 


-5.97 TO -2 


-1.54 TO -2 


dl 


0.071 


0.0 


a 2 


1.046 


0.0 


a 3 


4.97 TO" 3 


0.0 


04 


-4.38 TO -2 


-1.37T0- 2 


a 5 


-4.90 TO- 3 


0.0 


a 6 


-5.47 T0~ 2 


-1.37 TO" 2 



Table 2: The values of the effective wilson coefficient Cf and a, with the \x = mi, = 5.0 GeV, 
m c = 1.45 GeV in NDR scheme at NLO calculation. The coefficients a 2 i and a2i-i are defined 
by a 2i -i = C e iU + C 2 { f /N c , a 2i = C 2 { f + C e 2 {[ 1 /N c and we have taken N c = 3. 
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(pole/pole) at Tree level 
(pole/const.) at Tree level 



(pole/pole) with Penguin 
(pole/const.) with Penguin 



Browder et al. 15 



Hwang and Kim 1 [18 



Hwang and Kim 2 flT9"| 



Capstick and Godfrey p0| 
Dominguez plfl 
UKQCD || 
BLS [El 



MILC m 



WA75 [25 



CLEO 1 |26 



CLEO 2 pi 



BES |28 



E653 M 



f D . (MeV) 



322 ± 76 
322 ± 76 



336 ± 79 
336 ± 79 



243 ± 70 



362 ± 15 
336 ± 13 



fn. (MeV) 



196 ±41 
214 ±45 



219 ±46 
239 ± 50 



277 ± 77 



309 ± 15 
399 ± 24 



290 ± 20 
222 ± 48 
212±|±f 
230 ± 7 ± 35 
199 ± 8±t?±o° 



238 ± 47 ± 21 ± 48 
282 ± 30 ± 43 ± 34 
280 ± 19 ± 28 ± 34 

430^ ± 40 
190 ± 34 ± 20 ± 26 



fD* I fD 3 



1.64 ±0.48 
1.51 ±0.44 



1.53 ±0.45 
1.41 ±0.41 



0.88 ±0.35 



1.17±0.02 
0.84 ±0.03 



Table 3: The obtained values of fjj* (MeV) and fo 3 (MeV), and their ratio /oj//d s , and the 
results from other theoretical calculations and existing experimental results. Here we refered 
the corrected fn s values [34| for the experimental data [^5) - [EJJ . 
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